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ABSTRACT 


By a re-definition of the standard coefficients of form, 
wherein displacement is weighted by its distribution in depth, 
the analysis of wave making resistance is simplified. Using 
the Taylor Standard Series as a data base, regression analyses 
are employed to develop a universal wave drag function of only 
one variable for each speed length ratio. 


Form drag is decomposed into components which facilitate 
analysis. Prediction formulae are developed graphically from 
the Taylor Standard Series data. 


Using the prediction formulae developed, wave drag and 
form drag are predicted for a number of sample ships and are 
compared to that predicted by Gertler's reanalysis of the 


Taylor Standard Series. Similar comparisons are made for a 
half dozen models which are unlike the Taylor Standard Series 
parent with their individual towing tank test results. Dif- 


ferences and possible sources of error are discussed. 
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NOMENCLATURE 


sectional area 
sectional area weighted by its distribution in depth 
midship section area 


midship section area weighted by its distribution in 
depth 


free wave amplitude spectrum 
beam at midship section 

beam to draft ratio 
block coefficient = V/LBT 

form drag coefficient = R,/1/2 ov7s 


skin friction error coefficient = Ro/1/2 pv*L* 


coefficient of the form drag coefficient proportional 
to displacement = R,/1/2 pv*L* 


midship section coefficient = A_/BT 


midship section coefficient of the weighted sectional 
area = Ax /BT 


longitudinal prismatic coefficient = V/A, L 
weighted longitudinal prismatic coefficient = V*/A*L 
midshiv section coefficient = A,./BT 

residuary resistance coefficient = Rp/1/2 pv*s 
residuary resistance coefficient at V/YL = .5 

wetted surface coefficient = S//VL 

Pe anercic Goerficient = V/L’ 


weighted volumetric coefficient = V*/L° 


wave drag coefficient = R/1/2 pv*L* 
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displaced volume 

displaced volume weighted by its distribution in depth 
boundary layer thickness 

Froude number = V//’gL, where V is in feet per second 
skin friction error function 

form drag slope function 

interim wave drag function = (F4) (1000) *? 

universal wave drag function 

gravitational constant = 32.2 ft/sec’ 

wave number = 21/i = g/V? 

length on the waterline 

longitudinal center of buoyancy 


longitudinal center of the weighted displaced volume, 
Vx 


length to draft ratio 

wave length 

velocity potential of a 2-D wave 

velocity potential of a 2-D wave at z=0 

form drag 

residuary resistance 

wave drag 

skin friction resistance error 

form drag component proportional to displacement 
wetted surface 


draft at the midship section 
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w(x, Z) 


Saar 


direction of propagation relative to ship heading 
speed length ratio (L in ft., V in kts.) 

ship speed (in knots or ft. per second) 

offset of hull form 


half width of separated wake 
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I. INTRODUCTION 


Although many phenomena encountered in Naval Architecture 
have been reduced to relatively simple prediction formulae 
with reasonable engineering accuracy, there are as many which 
do not admit to any simple analysis. Among the most notorious 
are those pheromena which are responsible for the residuary 
resistance of ships. Withstanding centuries of scrutiny and 
data taking, they have steadfastly refused to collapse into 
any reasonably simple relationships involving the standard 
coefficients of ship hull form. 

During the 18th century, Chapman [3] built ship model 
towing tanks and developed scaling laws aii synthesis 
models from his observations. Unable to sort out all the 
phenomena involved in ship resistance with his ship model ex- 
periments and other observations, he complained, "For ships, 
we have to fear an infinity of bad qualities of the greatest 
consequence, which we are never sure Of being able to remove, 
without understanding the theory." 

During the last century since the publication of Froude's 
paper [13] in 1868, ship model basins have been constructed in 
many parts of the world for the purpose of measuring the re- 
Siduary component of ship resistance. 

Although we have probably sorted out most of the phenomena 
involved in ship resistance and have explained them with the- 


Ories, we have not yet succeeded in predicting residuary drag 
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in terms of standard ship form coefficients. This is due not 
only to the mathematical complexity of the theory, but also to 
the quantities required to be measured by the theory which are 
in general not those measured by standard coefficients of ship 
hull form. The difficulty is best expressed by Weinblum [26]. 
"... the wave resistance depends, to a first approximation, 
upon a complicated function of the surface slope in the longi- 
tudinal direction, i.e., on derivatives. On the other hand, 
the most commonly used hull coefficients are integrals which, 
even when kept. constant, still admit of very wide variations 
of the slopes. We realize now why the solution of the basic 
problem of the model basSins mentioned above - to establish the 
resistance as a function of the form - remains almost hopeless 
as long as the Ship surfaces (or at least their most important 
features) are not defined in a rigorous way by mathematical 
expressions." 

It is the purpose of this thesis to do just that, i.e., 
to relate residuary drag to hull form. The success of this 
endeavor depends on two considerations. The first is that the 
description of the hull form adopted for the thesis contains 
parameters which measure the distribution of displacement in 
depth and its effect on the velocity potential with Froude 
number, in place of some of the standard coefficients of form. 


The second is the word "almost" in Weinblum's quote. 
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IIT. WAVE DRAG 


The wave drag expressed in terms of the generated free 


wave amplitude spectrum, A(0), will be: 


+7/2 
R, =: 5 TeV? | [A(6)]*cos*6dé (1) 
-T/2 


The cos*@ term represents the weighting of the contribution of 
each part of the wave Spectrum to wave drag with the relative 
direction of propagation. The cos*@ weighting function indi- 
cates that the major contribution to wave drag will be from 
those waves travelling in the same direction as the ship (0=0). 
A Kelvin wave pattern is illustrated in Figure 1 (taken from 
reference 22) which shows ship generated waves propagating in 
various directions. In this illustration, it is the transverse 
waves following the ship which radiate most of the energy away 
from the ship. In view of the relatively small contribution 
to wave drag from the shorter diverging waves, it 1S a reason- 
able simplification to represent the wave system causing wave 
drag by the wave length and wave number of the predominate 
transverse wave. 

Several analytical methods have been derived predicting 
the wave drag of slender bodies. Mitchell's [14] integral is 
probably the method most often cited. Rather than standard 
coefficients of form, Mitchell's integral measures the form of 


the hull in terms of longitudinal slopes of the surface of the 
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hull at each point on the hull. The pressure change due to 


the resulting wave system is integrated over the ship to de- 


termine the total drag force. The coordinate system is shown 


in Figure 2. 


Using Mitchell's integral, wave drag may be expressed as 


(taken from reference 26): 


R = -2) | 6p ay dz = -2) [ap (32) axaz 


ff 96, ,3 
vad | pel (52) axdz 





or 
40g * 2 2 
R = =ES— | (1? + J?) A 
iE 
where 
+e K(x) _A*gz 
I = Eos ¥ cos (49%) axaz 
OX V 
ae 
and 
+e K(x) _ A*gz 
_ oy 2 ax 
J | (as) e sin(—3 )dxdz 
en 0 


where z = K(x) is the profile of the fairbody. 


(2) 


3) 


It may be seen from the foregoing that Mitchell's integral 
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integrates the slope of a given waterline, oy, weighted by a 


maACTOL, 


A?gz Dee 
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e | a (4) 


where F = Sy is a local Froude number based on the ship speed 
and the depth of the waterline, over the length of the ship. 
The area of each waterline is then weighted by its speed and 
depth below the free surface. The weighted waterline areas 
are in turn integrated over the draft, which weights each in- 
crement of disvlacement by its depth. 

The "effective" displacement of a ship which contributes 
to the wave making velocity potential is therefore something 
less than the real displacement by the weighting factor in 
Metehell's integral. If the lines of two vessels of equal 
real displacement are similar except for proportional changes 
in beam to draft ratio, the smaller beam to draft ratio should 
have less effective displacement and less wave drag than the 
vessel with the larger beam to draft ratio. The difference in 
the effect of the displacement at various depths is illustrated 
Poe Figure 3. 

In addition, due to the wave number, K, in the exponent, 
it is expected that the magnitude of any effect on wave drag 
due to the distribution of displacement in depth will also 


depend on speed length ratio or Froude number. For high speeds 
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So 


where the wave length of the ship generated waves is approach- 
ing the waterline length of the ship, the change in velocity 
potential over the draft of the ship will become relatively 
small for ships with large length to draft ratios. 

The variation of the velocity potential with depth may be 


written as: 


$(z) = o(z=0)e "7 = ge (5) 
where the wave number, K = a = <n and A is the wavelength of 


the transverse waves generated by the ship. 
At high speed length ratios, 


X approaches L, and 


K approaches af, 


The speed length ratio, 


V (knoxs) V (knots) 
YL Jr 


approaches 


V 


~ [Seem Soe 


and the Froude number, F 


= Aes approaches ae 
/GL /gh 
V ui 
= 7 = es = .399 
Vg (21V 77g) f/2T 


In addition, the effect on the velocity potential, 
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oe approaches >5e ; 





At this high speed, the effect on the velocity potential 


at various values of z is: 





1 -1/2 
b(2=F) = 6,677? = 12086, (6) 
i et /E 06S 
mr5,068 ~ %o° oS (7) 
=D 
b (z= = $e ie 8826, (8) 


The speed of a displacement vessel corresponding to a 
Froude number of om is often referred to as the "hull speed" 
for that vessel. It may be seen from the foregoing calculations 
that the magni.tude of the velocity potential will decay in 
depth from the free surface at a rate which will reduce its 
value to one half that at the surface at a depth of about one 
ninth of the waterline length of the vessel at hull speed. 

As ship speed is decreased from a speed length ratio of 
1.34, i1.e.,-down from "hull speed", wave drag decreases sharp- 
ly. The wave number, however, increases rapidly as speed is 
decreased, thereby increasing the rate of decay of the expo- 
nential with depth. There exists, therefore, a speed range 
wherein the wave drag has not yet become negligibly small, but 


where the rate of decay with depth of the velocity potential 
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has become very strong. In this situation the effect of the 
distribution of displacement in depth on the wave drag may 
become very important. Figure 4 illustrates the change in the 
decay profile of the velocity potential with ship speed. 

The decay of the velocity potential with wave number and 
depth may be written: 

$(z) = ge 8? = $e 92/v" (9) 

To demonstrate the variation in the decay curve of the 
effect of displacement with depth, caused by length draft 
ratio, we may look at the amount by which the effective dis- 


placement is reduced at the keel. Taking Z=T, 


= 2 
~ gT/V 
O 


eo” (GL/V*) (T/L) 
O 


o(2) = ¢ = 9 


2 
ep et 1F*) (L/7) 


= (10) 


where F is the Froude number and L/T is the length to draft 
ratio. 

It may be seen from this approximation that a large length 
to draft ratio will flatten out the decay curve in a manner 
Similar to that of a large Froude number. Long shallow ships 
or high speeds will therefore result in the least influence of 


the distribution of displacement over depth on wave drag. For 
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either long shallow ships or high speeds, in other words, the 
effective displacement will approach the real displacement in 


its contribution to wave making. 
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Pi ee ORM DRAG 


Without viscosity, the water is allowed to slip past the 
hull of a ship, its interaction with the ship being limited to 
normal forces and wave radiation. D'Alembert's paradox allows 
the pressure distribution aft to produce the same normal forces 
as those acting against the forward part of the ship. With 
the introduction of even the smallest amount of viscosity, 
however, a no-Slip condition is imposed upon the skin of the 
ship. This would directly result in the shearing of viscous 
layers of fluid very close to the hull and the diffusion of 
vorticity outward at a rate proportional to _ the 
gradient of velocity at the hull. ieadalt lon, at 
high Reynold's numbers, vorticity is physically con- 
vected away from the hull by turbulent mixing. By these 
mechanisms, momentum is exchanged with the flow and the velo- 
city distribution is changed from that of the inviscid flow. 
For very high Reynold's numbers this region of high vorticity 
flux and momentum exchange is limited a a relatively thin 
boundary layer which grows from the bow to the stern and con- 
tinues aft as all or part of the wake. Figure 5 shows the 
boundary layer thickness and velocity profiles along the hull 
of a ship. The presence of the flow entrained in the boundary 
layer prevents the streamlines aft from closing as quickly as 
those ie ene inviscid case. As a result, the pressure is not 


able to increase aft back up to the value which would produce 
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normal forces aft with forward components to equal the after 
components of the normal forces on the forward part of the 
ship. The difference in longitudinal components of the normal 
forces resulting from a reduction in the pressure field aft by 
the boundary layer is the viscous pressure drag. 

If the curvature of the hull lines aft is too abrupt or 
if for any other reason the pressure gradient in the direction 
of flow is strongly adverse, the boundary layer flow may sepa- 
rate from the hull, leaving a region of reversed flow following 
the ship. The region of flow following the ship constitutes a 
momentum defect, which can be much larger than that due to the 
boundary layer itself. 

Very roughly, drag due to flow separation may be repre- 
sented in a two-dimensional velocity profile with the reference 


frame fixed on the ship, as shown in Figure 6. 


+00 +6 
D = [puiu,, - u)dy = | pu(v,, - u)dy (10a) 
—0€O -6 


From Figure 3, it can be seen that for a separated wake, the 
contribution to this integral from separation along is approx- 
imately 


[pv2ay. 


This quantity may be quite large. If flow separation 
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occurs, then the form drag will increase by this large amount 
and may increase to a value many times that of the form drag 
due only to an attached boundary layer. 

An additional component of drag which is related to the 
form of the hull is the difference in frictional resistance 
between the wetted surface of the ship in a plane uniform flow 
and that in the curved form of the hull. Because flow is 
forced around the ship in a curved path, the pressure distri- 
bution causes an increase in the velocity near midships and a 
decrease near the ends. The increase in velocity around the 
midship region has a greater effect on the frictional drag 
than the reduction in velocity at the ends. The result is an 
overall increase in frictional resistance over that of plane 
uniform flow. In addition to the velocity field changes, the 
path is longer along a curved hull than along a flat plate of 
the same area. These effects are both more pronounced for 
shapes with high volumetric coefficients than for very slender 
forms with low volumetric coefficients. Since this friction 
drag is pelaced to the form of the hull, it is included in 
form drag. 

If form drag were plotted as a function of volumetric co- 
efficient and then extrapolated to zero volumetric coefficient, 
the value of the drag at the resulting intercept point, al- 
though included in form drag by conventional test procedures, 


cannot be due to the amount of displacement contained in the 
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hull form because the hull can have no thickness at a zero 
value of volumetric coefficient. Since frictional resistance 
is not measured directly in the towing tank, but is calculated 
by whatever frictional drag formula is adopted by the test 
procedure, then the value of the form drag at zero volumetric 
coefficient constitutes primarily the difference between the 
value predicted by the adopted friction formula and that 
actually contributing to total drag. Where they are signifi- 
cant, any other errors in the measurement, calculation, and 


scaling of frictional drag would also end up in this category. 
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IV. THE INTERACTION BETWEEN 
WAVEMAKING AND FORM DRAG 


The thickness of the boundary layer is dependent on the 
pressure gradient along the path of the flow. The pressure in 
general increases up to the stem, decreases toward midships, 
increases up to the stern and decreases back to ambient in the 
wake. Flow is thereby decelerated ahead of the ship, acceler- 
ated from the stem to midships, decelerated from midships aft 
and accelerated slightly into the wake. The Eluid entrained 
in the boundary layer is then flowing against the pressure 
gradient from the midship section aft. The boundary layer 
will therefore thicken at a greater rate aft than it will in 
the accelerating flow forward. The pressure decrease from the 
stem to midships is greater than the pressure recovery from 
midships to the stern. This causes the viscous pressure drag 
mentioned in the previous section. It has the somewhat re- 
deeming virtue., however, that the adverse pressure gradient 
against which the fluid in the boundary layer must flow is 
decreased in a gentle and continuous manner from that which 
would exist with a complete recovery of the pressure at the 
eeern., 

All of the foregoing refers only to the pressure distri- 
bution which would exist without waves. When the pressure 
distribution of a wave system is superimposed upon the pressure 
distribution discussed in the previous paragraph, skepticism 


arises concerning the predictability of the interaction between 
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two completely different phenomena which are each in their own 
Peqa@terelatively unpredictable. 

One of the most insidious interactions lurking in the 
towing tank might be the initiation of flow separation by the 
pressure field of the wave system. If the flow were on the 
verge of separation near the stern of a model and if the model 
were then run at a speed corresponding to a hollow in the re- 
Sistance curve, the crest of a transverse wave would be posi- 
tioned at the stern, adding a pressure peak near the end of an 
already adverse pressure gradient, whereupon flow separation 
might occur. This sudden increase in drag might then fill in 
all or part of an anticipated hollow in the resistance curve. 
If the speed were then changed to that corresponding to a hump 
in the anticipated resistance curve, a transverse wave trough 
would be positioned near the stern, reducing the pressure and 
accelerating tne flow toward the stern. Any region of previ- 
ouSly separated flow might be partially or completely re- 
attached, suddenly reducing the drag and removing all or part 
of the expected hump in the drag curve. 

In addition to the effect of the wave system on flow 
separation, flow separation may in turn affect the wave system. 
Mr. P.N. Joubert, in his comments on reference 9, cites an 
example of what appeared to be a wave induced by flow separa- 
tion on the 12 meter yacht "Valiant". 


In this way, the resistance curve of a model test might 
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be drastically altered by an interaction effect. If the model 
were then scaled up to a full size ship, the very large Rey- 
nold's number might prevent flow separation even in the pres- 
ence of the same wave profile as was modelled. The resistance 
of the full size vessel might then be entirely different from 
its model test. 

Even without flow separation, Similar effects might be 
produced ' to a lesser degree simply by the influence of the 
wave system pressure field on the thickening and thinning of 
the boundary layer over the afterbody. R.L. Townsin [25] said 
that the reduction in the non-dimensional boundary layer thick- 
ness, 65/L, may be as great between an 8-foot model and a 28 
foot model as it is between the 28-foot model and the full 
size ship. It is also agreed in the introduction to reference 
25 that separation scale effects, if they exist, may appear 


over the range of model size. 
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V. THE DATA BASE 


The selection of a methodical series of ship model tests 
for use aS a data base for establishing trends and functional 
relationships between wave drag poeocner parameters is limited 
by several considerations. The first of these considerations 
is the absolute size of the model. 

Regardless of the causes, the average magnitude of the 
errors encountered in tank test results increases strongly 
with decreasing model size and Reynold's number. Figure 7 
shows this relationship. It is suspected by the author that 
irregularities in tank test data may result from the interac- 
tions between the wave system pressure field and the boundary 
layer along the afterbody, as discussed in the previous sec- 
tion. Uniformity in the behavior of the boundary layer is 
improved by increasing the Reynold's number and since the 
speed and kinematic viscosity are not available as control 
variables for this purpose, then increaSing the size of the 
model will produce the same result. To eoce that wave 
effects are measured ina flow regime similar to that of a 
full size ship, it is extremely important that any irregular- 
ities due to small model size be avoided. For this reason, 
all ship model test data for models less than 10 feet long 
were rejected as unreliable. 

When one requires the size of the model to be larger than 


10 feet, testing becomes extremely expensive and methodical 
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Wemerres £Or large models in recent years tend to include only 
small numbers of models. This brings us to the next consider- 
ation, that of statistical significance. Data populations 
with only a few variations in some parameters of interest may 
lead to results which, regardless of their correctness, will 
@eeestand Up under any statistical measure of credibility. An 
example of this problem is provided by the BSRA Trawler Series 
(Part I) [17], wherein a total of only nine models were built. 
Four models with varying slenderness at a constant beam to 
draft ratio were tested along with five models of varying beam 
to draft ratio at a constant slenderness. With no variation 
in prismatic coefficient and so few in other parameters, Mr. 
G.S. Selman commented, "It surely needs a lot of imagination 
to produce from four models such wonderful curves, as shown..." 
[18]. Unless criticism of this nature is to be invited, a 
larger data base is required. This methodical series could 
not be used, therefore, even with a model size of twelve and 

a half feet and a scale ratio of only one 1/12. 

When other parameters are desired, the BSRA Trawler Series 
becomes more attractive. Although the number of data points 
in any one parameter is limited, changes in LCB, C. and load 
condition are available if Part I and Part II [17 and 18] are 
both used. Valuable insights have been gained from multiple 
regression analysis of the entire series [17]. 


Series 60 [24] was investigated and even with large numbers 
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of large models, could not be used because the range of pris- 
matic coefficient did not include values sufficiently low to 
be applied to many vessels other than merchant ships. In 
addition, the parent lines of Series 60 possess shoulders 
which, although appropriate for merchant ships, might intro- 
duce the additional complexities of the fore and aft shoulder 
wave systems Similar to those of the double-wedge-shaped body 
investigated by Wigley [27]. 

After consideration of several other methodical series, 
the Taylor Standard Series [7 and 22] was selected for a number 
of reasons. The models are all about 20 feet long and average 
about a ton in displacement. The number of data points avail- 
able in the variation of volumetric coefficient, prismatic 
coefficient and speed length ratio are sufficiently large to 
clearly represent functional relationships. The lines of the 
parent model for the roRpiter Standard Series are presented in 
Figure 8. Both LCB location and midship section coefficient 
are held constant in the Taylor Standard Series. 

An additional attraction to the Taylor Series as a data 
base is the manner in which the sectional area curve and the 
lines of the parent are changed to produce the desired varia- 
tions in form coefficients. The Sectional area curve is 
defined by a fifth order polynomial which is related directly 
to the prismatic coefficient and which is proportioned for the 


variations in the series in a smooth reproducible manner. The 
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lines are changed from the parent graphically in a process 
consistent with the changes in the sectional area curve and 
produce very fair hull forms without shoulders or other areas 
Of high local curvature. 


Almost 200 models were tested from 1902 to 1914, of which 


120 were published in The Speed and Power of Ships [22]. Al- 


most 160 of these tests were later faired into the data pre- 


sented in the Gertler re-analysis of the original test data 
[7]. Although the data is old, it is not out of date for the 
purposes of this thesis. In fact, a single methodical series 
of the magnitude of the Taylor Standard Series will probably 
never again be attempted. 

The data published in the three editions of Taylor's The 
Speed and Power of Ships [22] included only beam to draft 
ratios of 2.25 (Series 21) and 3.75 (Series 22) between which 
interpolation was recommended. Series 20 contained the data 
for the intermediate value of beam to draft ratio of 2.92. 
Although Series 20 was completed prior to Series 21 and 22, it 
was, for some reason unknown to this author, not published 
along with Series 21 and 22:as an intermediate value from which 
to interpolate beam to draft ratio effects. In addition, Series 
20 was not extended in 1913 and 1914 as were Series 21 and 22. 
Taylor did not explain why this data from some 38 models was 
left out, but only mentions in passing that it confirmed to 


some extent the validity of linear interpolation between Series 
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ee and 22. 

"A Reanalysis of the Original Test Data for the Taylor 
Standard Series" [7] — published in 1954 by Morton Gertler. 
Series 20 was included in these data and was used to fair the 
data and provide an intermediate value of beam to draft ratio 
of 3.0. The original Taylor data was corrected for many ef- 
fects and updated for modern conventions in this reanalysis. 
For this reason, the data presented in Gertler's reanalysis 
were used in this thesis with only one reservation. The 
behavior of the curves of residuary drag coefficient with that 
value at V/YL = .5 subtracted out and then multiplied by the 
wetted surface plotted against volumetric coefficient for each 
beam to draft ratio, speed length ratio and prismatic coeffi- 
cient Beyveal a very unusual behavior for the curve representing 
the beam to draft ratio of 3.0. While the curves for the B/T 
of 2.25 and 3.75 behave similarly, the curve for the B/T of 
3.0 usually starts out shooting above the other two and then 
crossing over the other two at a higher value of volumetric 
coefficient. This strange behavior of the data for B/T of 3.0 
is evident at many values of prismatic coefficient and speed 
length ratio when plotted carefully on a large scale. Examples 
are plotted in mares 9, 10 and 11. Because of this unusual 
and unexplainable behavior, the data for B/T = 3.0 is not used 
in this thesis. It is possible that a difference in the be- 


havior between the data for Series 20 and that of Series 21 
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and 22 may reflect some important problem with the models or 
the tests that was suspected or recognized by Admiral Taylor 


himself and was for that reason omitted from his publications. 
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VI. WAVE DRAG IN THE DATA BASE 


Since the Gertler reanalysis of the Taylor Standard 
Series was chcsen as a data base, the range of length to draft 
ratio for which data was obtained will be examined to estimate 
the magnitude of the effect of the distribution of displace- 


ment in depth for reasonable wave numbers and speed length 








ratios. 
fee Cc C = ¥ (11) 
pm 
Ca =~ .923 for all Taylor Series models 
C 
BT V i V 
72. 7 3 = (12) 
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Cae 
L? = (-E) pr (13) 
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The variation in the velocity potential decay rate seen by 
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this data base is then 
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Since speed length ratios from 0.5 to 1.2 were included 
in the data base, the variation of the effects of the length 
to beam ratios included in the data base on the wave drag will 
be estimated. Froude number is proportional to the speed 


length ratio as follows: 
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At a VR/VL = 1.2, Froude number is .357. Then, 
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It may be seen that at this speed length ratio, there is 
very little variation over depth of the velocity potential for 
the larger length to draft ratio. For the smaller length to 
draft ratio, however, the velocity potential is almost reduced 


in half at the draft of the ship. A given redistribution in 
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depth of sectional area or displacement would then presumably 
cause a larger change in wave drag for the smaller length to 
beam ratio than for the larger. 


At a speed length ratio of 1.0, the Froude number is: 


r= (.2976)-K = .2976 23 
oe ® io AAG = ® ( ) 
Then, 
1 ye. _ -1/(.2976)7(49.2) _ 
ae) max! TE 1.0) = $e = .795o, (24) 
L Vv _ _ ~1/(.2976)7(12.9) _ 
o (15) nin! TE = 1.0) = $e = .3196, (25) 


Again it may be seen that the variation of $(z) is quite 
large for the smaller length to draft ratio. The variation of 
o(z) for the larger length to draft ratio, although greater at 
this speed, is still not a large fraction of o—- 


At a speed length ratio of .8, 
VK 
F = tre? 1 Oe: = .2381 (26) 
and 


if V_ _  .-1/(.2381)%(49.2) _ 
Ne max’ (8) = b6e = .6996. (27) 
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o( (A) ane Pani SCH, ae = -2556) (28) 
At this speed length ratio, the variation of ¢(z) is 
large even for the large length to draft ratio. At the same 
time, however, wave drag has decreased to a small value such 
that the absolute change in wave drag due to a given redistri- 


bution in displacement in depth would probably not be very 


large. 
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It is hypothesized that the residuary drag may be decom- 
posed into three components: wave drag, form drag and skin 


| friction error. 


Rp = Ro + R, ke (29) 
ae _ h 27 2 
where R = wave drag = =pV“‘L‘*C 

W 2 W 

C. = wave drag coefficient based on waterline length 
oa a oe. : 

Ry = form drag = 5 V~L Cy 

Cc = form drag coefficient based on length 

R_ = skin friction error = i ov?L2c 

O 2 O 

Cy = skin friction error coefficient based on length 


It is further hypothesized that the wave drag will be a 
power law function of the displacement, V*, of the ship weighted 
by its distribution in depth and an unknown function, Fo, of 
the longitudinal prismatic coefficient of the weighted dis- 


placement and the speed length ratio. 


= SV7L7C. | (30) 
vx nN * V 
C, = Ups) 1IF3(C,", 31 (22) 


where n = exponent to be determined 
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PF. = unknown function to be determined 
Because it was expected that n would be approximately 2.0, 


a preliminary function was defined such that: 


e Wii, 2 * 6 6YV 
= [ (=) }1F2*, ; ar) } (32) 


Fo could then be more easily determined and PF, could be related 


back to Fo when n became known more precisely. 
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In addition, it is expected that the form drag will be 
proportional to the real displacement, V, and also proportional 
to some unknown function, F,,of the beam to draft ratio and the 
real (unweighted) longitudinal prismatic coefficient. 
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The skin friction error is defined as the value of the 
form drag when it is extrapolated to zero volumetric coeffi- 
cient. It is defined in such a way that it represents that 
Mart Of the form drag which cannot depend on the volumetric 
coefficient of the hull form. It would seem at first glance 
that since the hull form can have no volume to shape, the skin 
friction error should not depend on any form coefficient what- 
eer. The other form coefficients are not, however, forced 
By the definition of skin friction error to be zero and it is 
expected that some of the errors included within the definition 
may vary with beam to draft ratio and prismatic coefficient. 
In addition, since the real skin friction is a function of the 
hull form, the skin friction error will in general also be a 
function of hull form as follows: 
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VIII. DEFINITION OF NEW VARIABLES 


Because the variation of the velocity potential of the 
wave system over depth depends upon both the Froude number and 
the depth itself, new variables must be defined which relate 
this dependency to the speed and form of the ship. In order 
to include the effect of the vertical distribution of dis- 
placement on the velocity potential with ship speed, a new 
displacement parameter is defined which is an integral over 
draft with the integrand weighted by depth and Froude number. 


The area of each section may be weighted as follows: 
* + = ~ 
a) = | y(z,x)e gz/V dz (39) 


The weighted displacement is then: 


X=L 

7S | A(x) ax (40) 
x=0 
Xx=L zZ=T 

Ve | [ | a eae “dzJdx (41) 
x=0 z=0 


where y(z,x) are the offsets determined graphically from ref- 


Givi z ; ; 
represents the weighting factor for the 


erence 7, and e. 
contribution of the local displaced volume to the disturbance 
of the velocity potential. 


It may be seen that, in general, 
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and that for high speeds, V”* approaches the real displacement. 
As discussed in the section on wave drag and the data base, 
all of the displaced volume contributes approximately the same 
toward the disturbance of the free surface at very high speeds. 

In a manner Similar to the derivation of standard coeffi- 
cients of form, a new set of weighted coefficients may be 
derived which include in them their relative contributions to 
the disturbance of the velocity potential. 


Define a weighted midship section area, 
x 1 = “ 
A = | y(z)e as dz = atx) at x = §% (43) 


Then a weighted longitudinal prismatic coefficient may be 





defined: 
= ve 
oe — ART (44) 
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The longitudinal center of weighted buoyancy is: 
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Although a variation in LCB was not desired and LCB is 
fixed amidships in the Taylor Standard Series, LCB* will, in 
general, change with speed. The Taylor Series ships tend to 
have U-shaped sections forward and V-shaped sections aft. At 
low speeds, more V* will be lost forward than aft and will 
result in a migration of LCB* aft at low speeds. To insure 
that this variation remained small enough to ignore, LCB* was 


defined and calculated. 
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IX. DATA REDUCTION 


Within the data base of the Taylor Standard Series, a 
representative subset was chosen with enough data points to 
falearly establish trends between parameters. The subset of 
ships contains 70 ships with the following range of variation 


of coefficients of form: 


Seem OOl, .002, .903, .004, .005, .006, .007 


Fractional offsets were determined graphically for each vari- 
ation of os EO@mEO stations and the 0, :1, .2, .3, -4, «5, .6, 
-8 and 1.0 waterlines. 

For those residuary drag curves in reference 7 which are 
Flat at the speed length ratio of 5, the residuary drag co- 
efficient at this speed, Cel-5), is assumed to be essentially 
form drag. The form drag coefficient, Cy, based on length 
Squared, may then be written as: 


z Ry 
i 1/2 0V7L- ae 


In reference 7, the wetted surface coefficient is defined as: 
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f° VS ee 
and 
mV 
C., = 73 (49) 


so that the combination of coefficients 


 , Salo3) S v. _ Rp(-5)- 
C,(.5)C.vC, = (173 pvzs) (TR) | [3 = Wrens 7-1 (50) 


is equal to the form drag coefficient, Ci: 
The form drag coefficient, Cys was then plotted against 
=, for various values of a and B/T. The form drag coefficient 


was found to be essentially linear in Cc. and was rewritten in 


the form 
@: (51) 


where Fy is simply the slope of the curve of form drag coeffi- 


cient as a function of Co. Fy is, however, an unknown function 
of CS and B/T. The form drag is then 
1 


ee B 
Ry = ZV IC CUP) (Cys 7) (52) 


The intercept of the curve of form dray coefficient versus 
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Cc. extrapolated to C=0 1s taken as the skin friction error 
coefficient based on length. So, 


| 
a (53) 
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where Fy is not a function of Cie but 1S an unknown function 


of oD and B/T. 


The skin friction error is then, 


(54) 


ri Go 


ae = ZeV"L? FOC, 
Examples of the plots of form drag coefficient versus Cc, are 
provided in Figures 12, 13 and 14. Examples of plots of the 
form drag coefficient versus Ss are provided in Figures 15, 
16 and 17. Although these plots include data for the beam to 
draft ratio of 3.0, these data were not used in the derivation 
of form drag functions. 

Linear regression analyses were done graphically to de- 
termine reasonable fits for the functions By gp and 


op! B/T). The following relationships were chosen: 
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Plots of Py versus Cy are provided for three values of 
fein Figure 18. 
Plots of Fa versus oc are also provided for three values 


of B/T in Figure 19 with some data points included. 
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A computer program was written to do the necessary calcu- 
lations and is presented in one version in Appendix A. Offsets 
were read and interpolated to provide offsets at .7 and .9 
waterlines which were not available in reference 7. Complete 
full size offsets were generated for ten stations and ten 
waterlines for each of the 70 ships. From these offsets, V, 
SF Coe x Cn and Co were calculated. All integrals were done 
using Simpson's rule. Calculations were done for vessels of 
190 feet in length. 

For each of the 70 ships the weighted form coefficients 
were calculated for 15 speeds, providing 1050 weighted forms 
from the 1050 ship-speed combinations. Although this thesis 
contains formulae for predicting form drag, they were pro- 
grammed only as a means of comparing predicted form drag to 
Cp (V/VYL = .5). To prevent the validity of one hypothesis 
| from depending on the validity of another hypothesis, the 
values of CR (V/7L =; .5) punched from the Taylor Series were 
used as form drag in the calculation of wave drag. F. was 


then calculated as 


= z S ae 
mec, - Cp (-5) 1 (>) (c_) | (57) 
The product of the first computer program was both printed 
and punched for further use. Printed outputs for a number of 


example ships are included in Appendix €.. For each of the 
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1050 ship-speed combinations, the following vital parameters 


were punched on cards for later analysis: 


| 


| 


NAME 





SLR 
CER 
CR5 
CV OL 
CVT 
CPS1l 
CPN 
SL2 
BDR 


BZ 


QUANTITY 
V/VL 
CR(V/VYL) 
hs) 
Cx x 10? 
C,, re IMO 
CPp* 

CP 

SV Als 

B/T 


Fo 


The printed output contains additional information, including 


Bo! Py 


S, V*, Froude number, LCB*, Cn and CX. 


B 





Beh 


X. THE REGRESSION ANALYSIS 


The Statistical Package for the Social Sciences [2] was 


used to conduct all regression analyses and correlation studies. 
This very powerful analytical tool provides extremely easy 
interfacing and an output containing every imaginable measure 
of statistical Significance. Simple regression, multiple re- 
gression and stepwise inclusion multiple regression processes 
are available. 

The 1050 ship-speed Pon marcrons were filtered to elimi- 
nate misleading data and shakey extrapolations. A ship with 


Cy, = .007 and C, = .7, for example, at a V/YL = 1.2 will have 


Pp 
a residuary drag coefficient so large that it would plot off 
the page in reference 7 if the extrapolation were attempted. 
Over 100 ship-speed combinations, although extrapolated as 
data in reference 6, were eliminated for these reasons. 


Correlation studies were conducted in an attempt to deter- 


mine n in the equation 


2 2) xy 2 | 
C., F4 (Cr) F. (Cr) (58) 
for the entire 1050 ship-~speed combinations. After investiga- 
ting various values of n by computer and calculating best values 
of n by hand, n was fixed at a value of 1.8 for all speeds, so 


that: 





EO 


== Eee) 7. (59) 
‘It should be mentioned at this point that a value of n = 1.8 
lvery slightly favors high speeds from V/VL = 1.0 to 1.2. For 
lower speeds, a somewhat lower value of n would have provided 
a better fit. For the purposes of this thesis, however, a 
Wenstant value of n = 1.8 was used for all speeds. 
Due to the irregular manner. with which the relationships 

between the form coefficients and residuary drag change with 
speed length ratio, a separate multiple regression analysis 
was run for each of the 14 speed length ratios from .55 to 1.2. 
In addition to the variable, CP*, a series of five sine func- 
tions of CP* was defined and made available for stepwise 
inclusion into each multiple regression analysis as necessary 
to obtain the desired fit. Since the values of F. were punched 


in the data rather than F3, it was necessary to compute F., as 


3 
a variable within the regression analysis program for every 
sample along with the terms in the Fourier sine series. Ap- 
pendix B.1 contains selected listings and output printouts for 
One of the 14 regression analyses runs. Variables were defined 
as follows: 


ES | CPx 


evsT =ecv* x 10° 


CPSV 


(CPS1 - .48)/.27 
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CPSA = sin(m1CPSV) 
CPSB = sin(27CPSV) 
CPSC = sin(3m7CPSV) 
CPSD = sin(4nCPSV) 
CPSE = sin(571CPSV) 

F, = F,(cvst) ** 


The Fourier sine series was limited to five terms to 

avoid the representation of data scatter in the final function. 
Each stepwise inclusion multiple regression analysis then re- 
gressed Fy with CPS1 and the five sine terms. A constant is 
automatically added into the regression as the zeroth order 


term. The final form of the universal function, F., is then: 


Bee CK + C 


3 CPsie+ CACPSA pee eros +. -CPSC 


it B Si 


+ CpCPSD f- CRCPSE (60) 
The resulting values of the constants are listed in Ap- 
pendix B.2 along with various measures of confidence and 
statistical fit. An example of a typical listing of SPSS con- 
htrol cards used is also provided in Appendix B.1 along with 
parts of the SPSS output for that particular run. It should 
be noted that the F, in the regression analyses computer pro- 


3 


gram is larger than the F, defined in previous sections by a 


5 





ey a 







factor of alooo) -- = 3.98 due to the definition of CVST, 


evsT = cv* x 10° (61) 


To use this universal function for the prediction of wave drag, 
F. must be divided by (1000) °? Om: 
ae 


mee ev") ° = F, (cvsT):?/(1000) -7 (62) 


Tables of F = F, (cv*) °° are provided in Appendix B.3. 
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. XI. WAVE DRAG CONCLUSIONS 

For each value of V/VL, F, constitutes an estimate of a 
universal wave drag function of Cs alone and represents wave 
i coefficient per unit weighted volumetric coefficient 

raised to the 1.8 power. With the use of the weighted coef- 
mectents of form, cy and car and with the implementation of 
the universal wave function, Fo, the representation of wave 
drag is simplified considerably. The entire Taylor Standard 
| Series, for example, is collapsed onto one sheet of paper as 
'represented in Figure 20. A family of curves of F. versus oS 
is plotted in this illustration with one curve for each value 
of V/VL. It is reassuring that the curves of F, versus ck for 
various values of V/YL show minima similar to those found 
| using the standard coefficients of form. Figure 21 illustrates 
these minima for the Taylor Standard Series data in reference 
. 

The value of the exponent of displacement (or C.,) would 
be expected to be 2.0 from Mitchell's [14] linear ship theory. 
Weinblum [26] states that wave drag waS approximately propor- 
tional to displacement to the 1.6 power. With the form drag 
subtracted out of residuary drag, it is concluded from the 
results of this thesis that wave drag 1S approximately propor- 
tional to the displacement raised to the 1.8 power. Thus half 


of the discrepancy between the theoretical value of 2.0 and 


the previously observed value of 1.6 has been explained and 
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removed. 

| Figure 22 illustrates the effect of both oe and of oe on 
| 
the residuary resistance. Some sample curves of wave drag 


| 
coefficient predicted by the universal wave drag function and 


| 
by the Taylor Standard Series are provided for comparison in 
end 1x C. From these curves it can be seen that wave drag 
ie predicted fairly well by the universal wave drag function 
in most cases. In almost every case, the humps and hollows 
m- reflected in the predicted wave drag curve. The earlier 
rise in the wave drag curve for high values of os and the 

‘later but steeper rise in the wave drag curve for low values 
of C_ are both well predicted by the universal wave drag 
function. 

| The wave drag of vessels with lines differing greatly 
from the Taylor Standard Series is not predicted as well as 

the wave drag of vessels in the Taylor Series. Some sample 
curves of predicted wave drag and wave drag from model tests 
are provided for comparison in Appendix D. In most cases, 

wave drag is overpredicted for values of V/VYL above 1.0 and 
underpredicted for values of V/VL below 1.0. There are several 
possible reasons for these errors. 
At high speeds, vessels with wide sections, full water- 

lines and flat buttock lines aft will derive more dynamic lift 


than vessels with deep sections, hollow waterlines and full 


buttock lines aft. As a result, the former will tend to squat 
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more than the latter at high speed. The vessel with the deep 
sections and full buttocks aft may have less waelatevey anag ac 
low speed than the vessel with wide sections and flat buttocks 
but will have more drag at high speed due to the additional 
wave drag caused by squat. 

The Taylor Standard Series parent has deep V-shaped sec- 
tions and full buttocks aft. The survey launch and the 5.5 
meter yacht "Antiope", on the other hand, have very straight 
buttocks and wide sections aft, providing more squatting sur- 
face than the other vessels. Although the destroyer and patrol 
boat models have V-shaped sections aft, they are not very deep 
Baa the buttocks are not badly curved. Prediction formulae 
based on the Taylor Series would therefore be expected to 
underpredict the resistance of the latter three vessels at low 
speed and to overpredict at high speed. It may be seen from 
Appendix D that the patrol boat is the only one of these four 
which does not. behave completely as expected. 

In addition to the differences in wave drag coefficient 
due to dynamic lift from the afterbody lines, other differences 
in the lines may also cause errors in the prediction of wave 
drag. 

Variations in model size and testing techniques may also 
explain some of the error in stadineellng the wave drag coeffi- 
cient of those models sampled which were not in the Taylor 


Series. As discussed previously, flow separation and wave 
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interactions may cause scale effects over a small range of 
model size. Reference 9 discusses scale effects in hydrody- 
namic testing in great detail. Many plots are presented of 
predicted resistance coefficient versus speed curves to pro- 
vide comparisons between larger and smaller geosims. It was 
noted that even for widely differing hull forms, from merchant 
ships to racing yachts, the smaller geosim has a larger re- 
Sistance coefficient at low speeds and the larger geosim has 

a larger resistance coefficient at higher speeds. The cross 
over point is a little above a V/YL of 1.0. This is shown 
roughly in Figure 23. If the wave drag coefficient prediction 
curves in Appendix D are compared to Figure 23, it will be 
seen that regardless of the causes, the errors show strong 
Similarities to model scale effects. The predicted wave drag 
coefficient curves, based on the regression analysis of the 
Taylor Standard Series, follow the general trend expected for 
the larger geosim. The model test resistance curves for most 
of the models having hull forms unlike the Taylor parent 
follow the general trend expected for the small geosim, even 
though they were actually 12.5 to 16 feet in length. 

Since the change in the drag curve predicted from the 
Taylor Series caused by afterbody lines which derive more dy- 
namic lift at high speed is somewhat similar to the change 
caused by testing a smaller model, i.e., higher drag at low 


Speed and lower drag at high speed, then it would be very 
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@eertacult to tell which effect, if any, was contributing to 
the underprediction at low speed and overprediction at high 
speed of any smaller models which have flat buttock lines. The 
survey launch may be in this category although the model was 
16 feet long, and it is more likely that the flat buttocks 
caused the wave drag to be overpredicted at high speeds. 

Since "Antiope" was 24.3 feet in waterline length and the 
destroyer was a 20 foot model, it can only be their relatively 
straight buttock lines which cause their wave drag coefficients 
to be underpredicted at low speed and overpredicted at high 
speed. 

Although it is difficult to tell from their lines, the 
trawlers may derive slightly more dynamic lift at high speed 
than the Taylor models which could explain the underprediction 
of their wave drag at low speed and overprediction at high 
Speed. In addition, the trawler models were 12.5 feet in 
length and may suffer some scale effect when compared to the 
20 foot Taylor models. 


Lines drawings are provided in Figures 8 and 25 through 29. 





XII. FORM DRAG CONCLUSIONS 


Plots of form drag coefficient, 


oe Vi 
jgpv?t? ~ “s°r‘7n ~ °°) YS <) 


as a function of C5 for three values of ce and three values of 
| B/T are provided with data points in Figures 15, 16 and 17 to 

- show how the data was faired. Similarly, plots of form drag 
Mee eFicient as a function of C., for three values of 5 and 
_three values of B/T are provided in Figures 12, 13 and 14. 
From these graphs, it is evident that form drag is very nearly 
proportional to both ce and C,. Although data for the B/T of 

| 3.0 is plotted, it was not used to derive the equations for F 
| and Fy for the reasons previously discussed in Chapter V. Form 
drag was therefore assumed to be linear in B/T between the 


| Values of 2.25 and 3.75. 


The resulting equations relating a and F 


1 to B/T and a 


are as follows: 


Fy = 0027 + (Cc, - MENT. 022. 4 . 02033 (5 Be 5] (64) 


rr 
lI 


14.0 + (= - 2.25) [4.0 + 20.7(C, - .5)] (65) 


The form drag coefficient is then 





BQ) = 


L?2 
C. = (C+ C,) (4) (66) 


Miner c Cy is defined by Equation (38) and C, is defined by 


| 
fguation (36), or 


a 


. L* 
| Ce = (FL + F\CL) CS) (67) 
With each of the plots comparing predicted wave drag to 
model test results is listed the predicted form drag and that 
of the model test. Figure 19 also contains data points which 


represent the value of F, calculated from the equation for Fo 


1 
}and the residuary resistance coefficient at V/VYL = .5 from 
reference 7. The plotting of these calculated "data" points 
for Fy 1s not a complete demonstration of the accuracy of the 
equation Since it assumes the accuracy of the equation for Fo: 
The best evaluation of the accuracy of the equations for Fy 
and Fy is probably the direct comparison of the predicted value 
to the value of the residuary resistance coefficient at a V/VL 
of .5 from reference 7. 
Form drag iS predicted reasonably well for all vessels in 
the Taylor Standard Series as may be seen from Appendix C. 
Prediction is best for those vessels with volumetric coeffi- 
Cients in the middle of the range of data, approximately .004. 


For heavier vessels, form drag is underpredicted increas- 


ingly with volumetric coefficient. This effect is simply the 
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manifestation of the slight nonlinearity of form drag with 
G@mepiacement. The effects of Oe and B/T seem to be adequately 
represented by the Equations (64), (65) and (66) for the ships 
lJof medium and heavy volumetric coefficients. 

For extremely fine vessels, the effects of os become 
slightly overpredicted. As a result, form drag for vessels 
with low and low Cc. is underpredicted but is overpredicted 
for vessels with high SS and the same Co. This is another 
manifestation of the slight nonlinearity of form drag with 
displacement. 

For hull forms differing from the Taylor lines, presented 
in Appendix D, it appears that form drag is predicted relatively 


well for slender hulls but not so well for vessels with higher 





‘volumetric coefficient. The old destroyer in Appendix D, for 
example, is predicted reasonably well with a volumetric coeffi- 
cient of .00139. The round bottom patrol boat in Appendix D is 
predicted less accurately with a volumetric coefficient of 


-00246. The survey launch in Appendix D is underpredicted by 


about 500% with a volumetric coefficient of .0043. Both trawler 
models are underpredicted by over 100% with volumetric coeffi- 
‘cients of .00691 and 00830. The form drag coefficient of the 
5.5 meter yacht "Antiope" with a volumetric coefficient of 
-00627 is predicted reasonably well as long as the fairbody 
beam to draft ratio is used. The form drag coefficient predicted 


by the computer program for "Antiope" was inaccurate because the 
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} 

full draft was automatically used to calculate the beam to 
draft ratio. A misleading beam to draft ratio of only 1.35 

| 

; , 

was then entered in the prediction formulae for Fy and Fos ig 


is interesting that "Antiope" is predicted so well when the 


other vessels of similar volumetric coefficient are so badly 


| 
underpredicted. 


Taylor Series models with high volumetric coefficient are 


‘also predicted much more accurately than the trawlers. Few 


‘commercial vessels whose feasibility depends on carrying capa- 


ety and acquisition poe can justify the fine lines of the 
‘Taylor Standard Series. Merchant ships, for example, almost 
IB ways have some parallel middle body and shoulders which may 
me nerare additional wave systems and induce flow separation. 
The trawlers have propeller apertures and may not have after- 
body lines as smooth as the Taylor models or the 5.5 meter 
yacht. The Taylor parent hull form is mathematical and care- 
Sally avoids harsh curvatures. It is difficult to conceive of 
any commercial hull form which would treat the boundary layer 
more gently than one with'‘the Taylor afterbody lines. A 5.5 
meter yacht is a racing machine, the design of which is not 
constrained by practical considerations such as propeller aper- 
tures and carrying capacity. It might be expected therefore 
that a racing yacht would have less form drag than a fishing 


trawler of similar volumetric coefficient. Lines drawings of 


these vessels are provided for comparison in Figures 25 through 
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1295 The heavier of the trawlers, it should be pointed out, is 


| 


‘completely outside the range of volumetric coefficient sampled 


- the regression analyses in this thesis. If the lines of 
the trawlers do not treat the boundary layer as gently as the 
5.5 meter yacht or the Taylor Series models, then it would take 
only a small region of separated flow to explain the 100% dif- 
ference in form drag coefficient between the trawlers and the 
other vessels. 

The 500% error in predicting the form drag of the survey 
| launch cannot be explained by subtle changes in the afterbody 
|} lines or by any inaccuracy in the prediction formula. With a 
volumetric coefficient of .0043, a reasonably slender vessel 
well within the range covered by the prediction formula, a 
Horm drag coefficient of 2.5 x Or 1s about five times too 
large and indicates that the boundary layer must have separated 
from the afterbody or else the transom must have been submerged 
deeper than the design waterline. 


Although the equations for Fy and F, are only artificial 


1 
formulas fitted to the data rather than based on physical laws, 


there are some interesting characteristics to the gross behav- 


1or of F_ andF.. 
Oo I. 
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.0027 + (C, - .5)[.021 + 02033 (z ~ 2.25)] 


ry 
il 


14.9 + (3 - 2.25) {4.0 + 20.7(C, - .5)1 
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For ships with low values of 7 Neate, Skil LrLcrion, error 
becomes a constant independent of beam to draft ratio. For 

| =. 

Ships with low values of beam to draft ratio near 2.25, on the 


‘other hand, Fy becomes a constant independent of prismatic co- 


efficient. 
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XIII. WAVE DRAG RECOMMENDATIONS 


The range of data used in this thesis is limited to that 





which was available in the time available. The 1050 ship-speed 
combinations generated, while a substantial number, included 
speed length ratios only through 1.2. It is recommended that 
this effort be extended through higher speeds. If higher 
Speeds are attempted using the Taylor Standard Series as a 
data base, it will be found that the Ore of ce for which 
data is available is narrowed considerably from that of the 
lower speeds, eliminating the fuller ships. It is likely, 
therefore, that any investigation of higher ship speeds will 
have to be done using either destroyer type hull forms such as 
Series 64 or a high speed Trawler Series. 

In addition to the limits on speed length ratio, the 1050 
ship-speed combinations generated for this thesis included, 
for most speeds, only those values of C., between .001 and .006. 
This range of Cc. is equivalent to a range of displacement 
length ratio of 28.6 through 171.4, i.e., very slender ships. 
These values cf displacement length ratio are limited, for the 
most part, to military vessels. Extending the range of dis- 
placement length ratio to 400 or C. to a value of .014 would 
include most interests in the field of naval architecture. 

Tf the sectional area curve is normalized to 1.0, then 
the displacement is non-dimensionalized and becomes equal to 


m=. if f= 1.0 and A. = 1.0, then: 
P x 
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x=L iL 
a eX) _ = 
Y= | A(x) ax | AL ox = CS = Ho (68 ) 
x=0 0 oe 


where UU is the zeroth moment of the sectional area curve. 
The longitudinal prismatic coefficient may be thought of as a 
zeroth moment of the sectional area curve. Sima larly> a first 


moment of the sectional area curve may be defined as: 
x=L x=1 


_ xA (x) _ = e 
| xA(x)dx = | RAG = (C,,) (LCB) mas ae a HU LCB (69) 


t 


} 


| =() x=0 


and LCB becomes the ratio of the first moment to the zeroth 


‘moment of the sectional area curve. The second moment will 


manen describe the fullness of the ends more specifically. The 


‘third moment will describe the asymmetry of the sectional area 


curve in greater detail. The second and third moments together 

may be used to describe the momental skewness of the sectional 
area curve. These and higher moments of the sectional area 
curve may be defined and used as higher order terms in regres- 

sion analyses. Weighted variables could be defined, as was 

done in this thesis, to eliminate the dependence on beam to 

draft ratio. With a sufficient number of terms, it would seem 

| that the hull form could be adequately described so that wave 
drag could be predicted even for vessels not geometrically 
Similar. 


In this thesis, only the gross first order effects of 


displacement and beam to draft ratio were collapsed out of the 
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data. Many higher order effects such as LCB were, by the use 
of the Taylor Series, carefully avoided. Although the wave 
drag of a hull should be predictable from a few well chosen 


gross parameters without addressing every detail of the hull, 


there are many variables evaluated in regression analyses of 


eT = = . 
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total drag {16 and 23] which change drastically in relative 
importance. Reference 1 discusses the relative contributions 
of these parameters in great detail and presents them all in 

a single compact illustration, re-drawn as Figure 24. From 
the behavior of the correlation of these parameters, it may be 
anticipated that the variables will have to be very carefully 
chosen in order to take into account differences in hull lines 


which may have the same basic form coefficients. 
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NOTATION FOR FIGURE 24 


total section area at the AP 

total section area at the FP 

after perpendicular 

total section area at the station of maximum area 
beam at the station of maximum section area on L 


WL 


V/ (Ly Ay) 


pee / AV“ 
AL / (BT) 
A,/A,, 
ey A. 


position of the longitudinal center of buoyancy measured 
from the FP 


forward perpendicular 


angle of the buttock at 1/4 T 


at station Ly. /20 forward 
of the AP 


W 


entrance angle measured between the centerline and a 
tangent to the designer's waterline at the FP 


ship or model length measured on the designer's waterline 
total resistance of the hull in pounds 

wetted surface 

s/y2/3 


draft measured from the designer's waterline to the lowest 
point of the station of maximum area 


ship speed in knots 
displacement in tons 
volume of displacement in cubic feet 


N/A 01T, )3, displacement-length ratio 


L 
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FIGURE 26.6. 
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. XIV. FORM DRAG RECOMMENDATIONS 
Due to time limitations and to the relatively small 


contribution of form drag to the residuary drag at high speeds, 


a computer regression analysis was not utilized to derive for- 


mulae for Fy and F,. The functional forms of Ee and Fy were 


it 


instead determined graphically and were limited to linear 
functions of os and B/T. Although these simple linear predic- 
tion formulae work surprisingly well, there are some data 


points plotted after the fact which are not well predicted. A 


regression analysis sampling all 315 residuary resistance co- 


efficients at a V/VL of .5 or less from reference 7 would 


surely pull in most of the wild points, particularly if higher 
order terms were made available for stepwise inclusion. 

If a regression analysis were set up to fit prediction 
formulae to form drag, it would take only a little additional 
programming to include other higher order effects such as LCB 


location in the analysis. 
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APPENDIX A 


Program Used for Calculation of Weighted 


Variables and F 


2 


for the 1050 Ship-Speed Combinations 


; 


cLOWZISA 


«> 


1189 107 boat Gato 


-Ul Sc2 302 -% BABA 


~ 





a Lo 








u 
C AN ANALYSIS OF RESICUARY ORAG 
— i ee 
1 OIMENSTON XSLOPCIL 11) sZ2ZSLOP(C 11411) ,0AC 11511) 
—~2 ~ ~~ OFMENSTON YYCLIs LIP eVYCL Le LL) pAACLTI LL) ,ACL1),AM{11),G0L1) ~ 
3 OIMENSTON CR(15,7) 
- 4 INTEGER VL,CVV,BC OO al or RI rr 
5 REAL M 
—on060€=~*~CrRK EAL LC 
7 REAL LBP 
~ 8 REAL LWL eee eerie. a 
9 INTEGER WL . 
UC GER SU 
11 BUG=.5 
See CLWL 19000 
13 LBP=190. 
14 —— -T#=4 2656 ee eee 
15 B=17.461 | 
sp 
16 READ(S,1) CCYY(WLsSTA) »WL=1,S),STA=1,11) 
—17———1- FORMAT (9F8. i ae aay os 
C 
“18 ~~~" READ(5,39) COCROVEs CVV) »VL=1515) sCVV=1y 7) ysis > aa 
19 39 FORMAT (15F5.2)_ 
eS - —— a a ee ee sk es _ 
20 DO 35 STA=1,11 
————-.-— 2 ete ae pao, SOS CRS 

















C ADOITION OF WATERL INES 9 AND ef TO THE TAYLGR CFFSETS 
sel ~~ 00 20 L=ly7? i a . _ 


22 WL=12-I 

Meee — VY(KL,STA)=YY(WL—-2,STA) ~~ 
24 20 CONTINUE 

ea Y(3,50)-YY(2.,StA} 
26 VVC 2¢STAD=(YY (1, STAD+YY(3_STA) 140.5 

—27- > HVS STADECYYO3 9 STADFYY(5 STA) ) 4005 nS 

C 

28 ~~. 35 CONTINUE a a co a Pe om er 
29 \ CVvv=l | 

aco 


. 


C 
—C PROPORTIONING OF FULL OFFSETS FRCM TAYLOR FRACTICNAL OFFSETS 
31 00 7 STA=1,11 





22 : OO 7 wl=1,11 vw: - SESS Se a as ae 
33 VOWEL» STAI=YY(WLE,STA)*B8*0.5 
—~34 —_——_ P- CONTINUE _— a a ee a 
C 


 —- €-> CALCULATION OF SLOPES AT EACH CFFSET REFERENCE PCINT 
3° OO 21 Wi=1l,11 





~~ 36 OO 21 STA=2,10 ~ aa ~~ 
37 21 XSLOPCWL »STAD=(Y Chl STAtL)- Yew STA 1))*05. O/LWL 

—~38°>———~-~~ ~DO 22 STAzlL,ll ~~~ aaa EEE 
39 DO 22 WL=251C 

~ 40 ~-22 ZSLOP(WLEeSTAI HC VYOWL Fl sSTAI-Y(WL-1,STAII*O5.0/T ~~~ an. yy 


41 DO 23 WL=loell 





Oe 


—— ee ee me ee ee a _o— +s - -_—_— — 


2 23 XSLOP(WL»1L)=2.0*XSLOP( WL »2)-XSLCP( WL, 3) 
3 DO 24 WL=1,Ll aed 
4 24 XSLCPC WL eo LL)=2.0*XSLOPE mL» 1CI-XSLOP CHL, 9) 
See DO 25 STAcI,11 © ~~ *° OO 














6 25 ZSLGPC1L»STADH2.08ZSLOP(2,STAI—-ZELOP(3e¢STA) 

LT DO 26 STA=1,11 sin io 2 

8 26 ZSLOP(LI»STAJ=H2.,08ZSLEOP(10,STAI-ZSLOCP(9,STA) 

Cc 

Yd VL=1 

| -_— a a See 

il BUGR=BUG+. 5 
C 

52 WRITE (6,19) BUGR 

3 19 FORMAT (81%, 'eeeeeeeseeeeceeeseseeeeessesee = 6SHIP NUMBER = ',F5.2, 

Cc? SSSSSESSSSESESSES OSSSHS SSS SSS SHSSOSSS SS seSe ses sess se se segesgess?) 

-———¢—-— ee es eo a” 

4 FR=z0.0 

5 GO TO 34 
Cc 

56 6 V=SL*SORT(L WL) *1.€89 

’ 

Cc 

-———€ —- CALCULATION OF EXPCAEANT FOR REIGKTING ENCREMENTAL AREAS §~— ee Te | oe a 

7 - FR232.2/ (V2) 
c . ae: 
C CALCULATION CF FROUCE NUMBER 

58 Fev/((32.-29LWi 094.5) * Sn eS a ee : ia 

Cc 

9 - - 34 CONTINUE : Pa Ohm © on mien rl oo eet om ae 

50 STA=l 

61 2 CONTINUE 
Cc 


C INTEGRATION OF WEIGEFTED INCREMENTAL AREAS OVER THE DRAFT FOR EACH SECTION 
62 OO 31 wl=l,1l 








S35 — 31 AACwWLeSTAD ZFYOCWL STA) FEXP(-FROCbL-1) *#T/10.0)%2.0 ” 
65 OUDC=0.0 
66 DO 32 wWi=2,10e2 : 
é? 32 EVNNZECVNNG AA (WL, STA) 
68 DO 33 wlL=3,9-2 
°69°— ~ ——33 OODO=O00D0t AA(WL se STA) i) a ah ia an Mess eae a a ee al 
70 AISTAID=A(T/3Z0-O) FC AAC L eSTAD +S .O08%EVNNG2.080COC+AA(LI,STA)) 
C 
C CALCULATION OF LCNGITUDINAL MCPEATS OF SECTIONAL AREAS 
Tl AM( STAI =A0 STA) ®FLCAT(STA) 
, c 
mT LE CER SGT oO oO GOT 
C 
C INTEGRATION OF INCREMENTAL GIRTHS OVER EACH HALF SECTION : ne at 
73 OO 27 WL=l,ll 
74 27 DAC kL» STAD =HASQKRT( LOF( XSLOP( Whe STA) I S826 (Z2SLCP (wh, STA) ) & 42) 
75 EVN=0.0 
oo - —— OOD0C=0.9 ie ee Mie) 1° 7... 0) =? oes 
77 DO 28 WL=2-10-2 
78 28 EVNZEVN*NA (WL, STA) “Hi aS 


19 DO 29 wWl2=3,9e2 
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80 29 ODOCZzO0ND+DAIWL STA) 


81 GISTAI=H(T/30-0)F(CA(L, STA) #4 OF EVN42.0* 0D000¢DA(IL,STA)) 
c 
MG CONTINU 0 
C 
83 STA=STA+tl 72 s me 
84 IF (STA.LE~L1) GO TC 2 
Cc 


C INTEGRATION OF SECTIGNAL AREAS OVER LENGTH FOR DISPLACEMENT 
MCIOO= 0.00 9 ee 








87 DO 3 STA=2,10-2 
88 3 EVEN=EVENtFA( STA) 
89 DO 4 STA=3,5992 
90 4 OOO=O0D+A( STA) 
— 9r DL SPL=(LBP/30.I*(ACLI4S. FEVERF2Z0900CHAC LI CO 
C 
Cc INTEGRATION OF MCMENTS GCF AREAS CVER LENGTH FGR LCB 
92 AT=C.0 
93 AMT=0.0 
94 00 S$ STA=1,11 
— 93 ~-— 9 AT=AT#+A(STA) ~ oo te 
96. 00 5 STA=1l,i,1 
97. ~ $§ AMT=AMT+AM(STA) : 
98 . LCBZCAMT/AT—6.)*°(LBP/LwL)*.1 
ae = oe eee = Ee z ee ate Sh RE se cy Oa es 
99 M=Lwl/(O1SPL**. 333333333) 
= =~ aC ee ne ee eee = = BAe meen ce ae PE ae a a ee 
C 
100 CP=DISPL/(A(6) FLL) 
101 CPS=CP 
Cc 
Lo2 Cv=DOISPL/(LWL**3) 
~ 103° -- ~—~ ~~ CvST=CV*1000.0 rr einen oe 
Cc 
10% CM=A(6)/(6*T) 
Gc a 
105 CB=CM*CP 
C ‘ 
~106—~———__ TF tFR.LE-0.02) GO TO 147 
Cc 
Cc F2 FOR CW BASED CN L¥¥2 | inns 
107 CCR=CR(UVL, CVV) 
108 CR5=CR(1.CVV) . an 
109 CRIVL»CVVI= O.OOLFCRI VE CVV) 
—..07-— F2=(CRIVEL, CVV) — CROEL,CVVI © St2 © COCELWE**ZI/(OITSPLII*e2) 0 ~~ 
Cc 
lil PUNCH 44 , St,CCR,CRS,CVST,CVT ,CPS,CPN, SL2Z2,80R,F2 
ll2 44 FORMAT (F422 2X 2F Gece 2X eGo le 2keF5e Se 2Xe Fe 3p 2X pFb eo Gy 2X oF bn 49 2X, 
CF6.4, 2X eF Ge 29 2Xe F622) 
C 
~113-~-—~-~——~ WRITE (6,23) SteFeOISPLsCP,CheCOrsCVeCCRetCBeFQ2 OO OO aaa 
ALS L3 FORMAT (90O* 6F 4.292% FS ode IX oF Ee Ce 3X 9 F904 oe 3K FD eG 3X yp FV 0% 3X FI 2G, 3K 


Ce F644 3X0F 904 SXF IZ. 6) 


tm 





See 


SHG = 








115 SL=SL+.05 _~ 
116 VL=VLel | = 
117 1F (SL.LE~1L.2) GO TO 6 
_——to6U™”™*~D™C™~« MmmCmhLT0)lm3860lUmUmUmUUUCU i re a eC ae ee ss 
C 
119 14 CONTINUE 
C 


C INTEGRATION OF GIRTFS QVER LENGTF FOR WETTED SUPFACE 
120 EVE=0.0 














Ta) 0 ES 0P5T 6 a a ce 
122 DO 16 STA=2510%2 
123 16 EVE=EVE*G( STA) 
124 DO 17 STA=3+902 
125 17 GO=CO+Gl( STA) 2 
126 $=0.06667*LBP9(G(1)44.%EVE*2.*CC*G(11)) 
(e CORRECTION TO WETTEC SURFACE FCR SKEG 
127 S$=S-T*Lwlh *(0.3372-0.438*CP) 
GC 
128 WS=S/SORT(DISPL *LWL) 
C 
ests NCP 
130... DISP=DISPL 
|) CVT=CV*1000.0 “a 
132 VR=CV 
133 “  -* §$L2=S/ (LW F*2) ay ee nC aicaiae °°. Cie ts 
134 BOR=B/T 
— 135 ~-— —_____ FO= .09027*+(CP-—.5) f -O2Z1#(30R—-2.253%.C020333)2 °° 7° Tee 
136 FO=FO*0.001 
137 Fl=14.¢ (BDR-2.25)4(4.4+20.7%(CP-.5)) 
138 Fl=F1*0.001 
139 FF=VR*F1l ¢ FO 
c 
— 140 — ~~  WRITECO,LLIILWE eBeTeDISPL»CP CM” ~~ ma ec 2 ee ar 
141 LL FORMATO ®—* g*lLwl =" F762 95X09 B = 8p F602, 5X9 TF = 8, FO02,5X,°OISPL = 
Cc *,FRE~1,5X,°CP = ©,F97S fe OX, SCH = *,F9.5) 
C a 
142 WRITE (6,12) CBeCVeSehS 
143 12 FORMATO ®—%,9CB = %9F Gee SXe* CW = 8p F965 SXe*S = 8 gFIOL0,5X,°CS = ° 
on.  —CrF9.4) ee ap — sax ee ee ny eae Sees AO Ste : 
C 
144 WRITE(6,37) SL2,FF 
145 37 FORMAT(C®°—-° ,*S/LWL9*2 = *,F10.8,10X, °FO + FILYCISPL/A(LWLe*3) =°,FLO. 
C8) 
C 
—"146 ~~ ~-> WRITE(6,8) Neate eee lf an is - 
147 8 FORMAT ®-*% g®SLRe op SX F® oSXe*CISEL 8° ,5X,"CP $°,8X,°CM £°,9X,°CB &? 
CoIXe*CV &°,5X-"* CCR *,3X,°LCE *° ,LOX,°F2Z") 
C 
148 GO TO 6 
C 
—~t1469 ~———36 CONTINUE ee eee a = ee 
150 Cvv=CVVel 


C DIMENSIONS OF ADODITICNAL SHIPS TC BE CALCULATED aa 
151 BUG=BUG+1.0 





152 
b>3 
154 
5D 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
L169 


170 
ivi 


172 
173 


100 


SENTRY 


T=6.5851 
B=24.6941 

IF (BUG.LE.2.0) 
T=8.065 
B=30.24%4 

IF (BUG.LE.3.0) 
T=9 2.3127 
B=34.9227 

IF ({(BUG.LE-4%.0) 
T=10.4119 
B=39.0447 

iF (BUG.LE.5.0) 
T=11.406 
B=42.7T 

IF (BUG.LE.6.0) 
T=12-3196 
6=4E.1984 

IF (BUG.LE.7.0) 


WRI TE(6,100) 
FORMAT(*1%,* END 


STOP 
END 


GC 


GC 


GG 


GC 


GC 


GC 


CF 


To 15 


TO 15 


TO 15 


VOnL5 


TO 15 


TO 15 


PROGRAM® ) 


=O D = 


| 
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APPENDIX B.1l 


Listing of Control Cards and Some Output Examples 


for One Regression Run 





1. 


I,.G AIGWAIIA 





smnod t 


bos abteS Lloxysgnod 3a 


mus Aole2sitpel en Jo? 


qweo— 
~ eee ————— 


= Ibi I= 


ixample of first page of printout for v/VL = 1.0 run 








PATISTICAL PACKAGE POR THF SOCIAL SCIENCES SPSSH - RPLBASE 6.02 03/25/76 

| eee _ eee 

P SPACE ALLOCATION POR THIS RUN.. 

7 TOTAL AMOUNT REQUESTED 80000 BYTES 

_ 

DEPAULT TRANSPACE ALLOCATION 10000 BYTES 

= i i a ae Se ee es eee eee ee eee a ae 
& : “7 MAX NO OF TRANSPORMATLONS PERMITTFD 100 

MAX NO OP RECODE VALUES u00 

ia MAX NO OP ARITHMN.OR LOG. OPERATIONS 800 





RESULTING WORKSPACE ALLOCATYON 7 "JOO00 BYTES  —<—™—— 


RUN NAME MULTIPLE REGR®SSION ANALYSIS OF THE TAYLOR SERIES 
PILE NAHE DATA POR PREDICTION OP WAVE DRAG 
Sees tt s)06™6™”™C€CUMARTABLE LIST.)OUSL&,CR,CRS,CVST,CVT,CPS1,CP,SL2,BDR,P2 - 
VAR LABFLS SLB SPEFC LENGTH RATIO / 
—_— °° CR TAYLOR RESIDUARY DRAG COEPPICIENT TIMES T0600 7 ‘ 


CRS CR AT A SLR OF 0.5 
—_ °° ~ , CVST WEIGHTED VOLUMETRIC COEBPPICIENT TIMES 1000 / 
CVT ACTUAL VOLUMETRIC COEFFICIENT TIMFS 1000 / 

—_ CPS! WEIGHTED PRISMATIC POLYNOMIAL TERM IOP ORDER OWE / 
; cP ACTUAL PRISMATIC COEFPICIENT / 
ri ~~ SL2  £RATIO OF WETT2D SURPACE TO LENGTA SQUARED FJ = 
BDR BEAM TO DRA?T RATIO 
z ae FP2 CV STAR SOQUARZD TIMES CW BASED ONL / 

INPYL PORMAT PIXED (FU. 2, 2%, FU. 2,20, 24.2, 2%, 8 5- 302% F523, 2%, F6-4, 2X 
= F629, 2%, F6e0, 2% ¢F 8.2, 2X, FS. 2) 





a i ee a = me 8 ge i ee * 


| ACCIRDING TO YOUR INPUT PORMAT, VARIABLES ARE TO BE READ AS POLLOWS 





| ___—sSWARIABLE PORSAT RECORD —-COLUANS i 
| 
) SLR P 4, 2 1 1- 64 
ee CRF 10 Se 
CRS P 4. 2 13- 16 
a. . cyst Pasir) 1°" 19-23 : oS tat 
CVT P 5. 3 1  26- 39 
- cPs1 P 6. 4 1 33- 38 
cP P 6. 4 1 1-46 
———— Pe (ie at 64-02 — 1-09-54 oS epee 2 
BDR F 4. 2 1 57- 60 
ae P2 P 6. 2 1 63-68 


HE INPOT FORMAT PROVIDES FOR 10 VARIABLES. 10 WILL BY RZAD 
T PROVIDES POR 1 RZCORDS ('CARNS') PER CASE. A MAKINGM OP 68 *COLUSNS* ARE JSED ON A BECORD. 


(er er a Ae AL ee na tye em apres er ed ee eee ee - 





COMNPUTP P3=P2* (CVST** (9. 2)) 
COMPUTE CPSV2(CPS1-.48) /. 27 
COMPUTE CPSA2SIN (3. 18*CPSV) 
COMPUTE CPSR=SIN (3. 18*CP3V%2.) 
CON PUTP CPSC=SIN(3. 18*CPSVE 3.) 
|. _ zzz _ . COMPUTE ss CPSD=SIN (3.14"CPS Ven.) Im : 
COMPUTE CPSZ=SIN(3. 18*CPSV*5.) 
8 OF CASES 66 
REGRESSION VARIABLES =F3,CPS1, 


CPSA,CPSB,CPSC,CPSD,CPS® / 
REGRESSION=P3 WITH CPS1(1), 
CPSA(1),CPSB(1) ,CPSC(1),CPSD(1) ,CPS&(1), 
RESID=0 / 

STATISTICS . ALL 


oN ee ee we re ee ee + tee 





See 


Example first page of printout for V/VL = 1.2 run 


STATISTICAL PACKAGE FCR THE SOCIAL SCIENCES SESSH - RELEASE 6.02 03/25/76 





SPACE ALLOCATICN FOR THIS RUN... 


TOTAL AMOUNT REQUESTED % — 80000 BYTES 





DESAULT TRANSPACE ALLOCATION 10000 AYTFPS 
— "MAX NC OP TRANSPCRMATIONS PERMIITEL 100 ‘ ee ee 
_— MAX NO CP RECODE VALITFES 400 
MAX NO OP ARITHM.OR LOG.OPERATIONS 800 a oa aes 





| RESULTING WORKSPACE ALLOCATICY 70000 BYTES 


— RUN NAME MULTIPLE REGRESSION ANALYSIS OP THE TAYLOR SERI®S 
BILE NAME DATA FOR PREDICTION OF WAVE DRAG 
—_— °°» VARTABLE LIST SLR,CR,CR® ,CVST,CVT,CPS1,CP, SL2, BDR, P2 Wee 
VAR LABELS SLR SPEED LENGTH RATIO / 


~~ CR — TAYLOR RESIDUARY DRAG CCEPFICIENT TYMPS TOO / 
CRS CSI AT A SLY OP 9.5 7 
eeetti(<i‘C:Ci‘i‘i‘i‘i‘i‘i‘i‘i‘i‘i CV ST)) 860 EIGHTED VCLIIMETRIC COEFFICIENT TIMES 1000 7 
CUT ACTUAL VOLUMETSIC COFPVZICIENT TINPS 1000 / 
EEE CPS1 WEIGHTED PRISMATIC FOLYNOMIAL TERM CF ORDER ONE / 
CP ACTUAL PRISMATIC COEPPICI®ENT 
TO SE DCWmRAATYSD OF WETTED SURPACE TO LENGTH SCUARED 7 
BDR FEAM TO DRAFT RATIO 
ar eal +’. CV STAR SOUAREL TIMPS CW BASFD ON L / 
: INPOT POPMAT PED CEG. 2 oe os 2 ote dees 2h, Pewee ekeF oe se lke Oe 2k, 


ere es | es 


P64, 2X, FEU, 2X, PU. 2, 2X, P74, 2) 





ACCORDING TO YOUR INPNT PCRMAT, VARIABLES ARE TO BE READ AS FOLLOWS 
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APPENDIX B.2 


Listing of the Resulting Regression Coefficients 


for Each of 14 Runs 


Coefficients are listed in column "B" 
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APPENDIX D 


Samples of Wave and Form Drag Prediction 


for Vessels Unlike the Taylor Standard Series 
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